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A new Bayesian image segmentation algorithm is proposed by combining a loopy belief 
propagation with an inverse real space renormalization group transformation to reduce 
the computational time. In results of our experiment, we observe that the proposed 
method can reduce the computational time to less than one-tenth of that taken by 
conventional Bayesian approaches. 
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Bayesian segmentation modeling based on Markov random helds (MRF’s) is one of 
the interesting research topics^^ Image segmentations are required to classify pixels in 
an observed image into several regions, and such segmentations are regarded as a kind 
of clustering of pixels. Markov random helds are regarded as classical spin systems in 
statistical mechanics^^ Loopy belief propagations (LBP’s) have been applied to con¬ 
struct certain practical algorithms for application in Bayesian image segmentations^’ 

In the present short note, we propose a new Bayesian image segmentation algorithm 
based on a RSRG transformation to reduce the computational time. 

We consider an image as dehned on a set of pixels arranged on a square grid graph 
(V, £). Here V^{i\i = 1,2, • • •, |V|} denotes the set of all the pixels and £ is the set of all 
the nearest-neighbour pairs of pixels {i,j}. The total numbers of elements in the sets V 
and £ are denoted by |V| and l^”!, respectively. The square grid graph has the periodic 
boundary conditions along the x- and ^-directions. The label at each pixel i is regarded 
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as a state variable, and it is denoted by aj. Each pixel i takes on the values of all the 
possible integers in the set Q={0,1, 2, • • •, g — 1} as its region label. The state vector 
of labels is represented by a = (aj|zGV) = (ai, 02 , • • •, a|v|)^- The prior probability of a 
labeling configuration a is assumed to be specified by a constant u as 

P(a) oc JJ exp(^^«ha,,„.), (1) 

up to the normalization constant. 

For the prior probability distribution P(a|u), we introduce the following RSRG 
transformation: 

exp(^^Q;^’'^hai,a3) « + ^a4,a3)) 

a2^QaA^Q 

(r = l,2,---,P), (2) 


where a^^'^=a. The square grid graph is defined as shown in Fig.1(b). 

Here VT) is the set of all the pixels and represents the set of all the nearest-neighbour 
pairs of pixels {i,j} in the graph Here, we remark that VT) is the subset of V. The 

prior probability distribution for a*^’')=(aj|zGVT)) on after 

the r-iteration of the RSRG transformation is expressed as 

n (3) 


up to the normalization constant. The transformation in Eq.(2) can be reduced to 
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The intensities of the red, green, and blue channels at pixel i in the observed 
image are regarded as state variables denoted by and df''^®, respec¬ 

tively. The random fields of red, green and blue intensities in the observed color im¬ 
age are then represented by the 3|V|-dimensional vector d={di,d 2 , ■ ■ , where 

di = df'’®®'^, df'"®)"*". The state variables df®*^, dp®®" and df''^® at each pixel i can 

take any real numbers in the interval (— 00 , -|-oo). It is to be noted that, in the above 
generative process of natural color images, we assign each pixel z to a labeling state Oj. 
For example, if the pixel i is in the labeling state Oj = ^ (C^Q), then its color intensity 
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Fig. 1. 
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Block spin transformation, (a) Graphical representation of Eq.(2). (b) 


(r = 1, 2, • • •) in the coarse-graining process. 


(y(r)^£-(r)) 


vector di is assumed to be generated from the following Gaussian distribution: 


^ V det(2xC(q) "^P( ■ - -"K)))' (5) 

In other words, the labeling state a* specifies the distribution within the set 
{ 5 f(dj|^)|^GQ} that generates a color intensity vector di. As mentioned in the pre¬ 
vious section, we introduce the labeling state variable o* at each pixel i as a Potts spin 
variable in statistical mechanics. 

After setting i? as a positive integer, we construct a 3|V*'^^|-dimensional vector 
= {di\iEV^^^) from our observed color image d. We remark that di in the coarse¬ 
grained image d^^^ is always the same as di in the observed color image d in our 
scheme. By using the Bayes formula, we introduce the posterior probability distribution 
for as follows: 

oc ( JJ ( Yl (6) 


up to the normalization constant. The estimates and C(^)|^gQ} of 

and C(^)|^gQ}. are approximately computed by means of the LBP algorithm 

for the conditional maximum entropy framework^b After obtaining these estimates, we 
calculate from by using the following inverse RSRG transformations of Eq. (4): 
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where a = 5;^°^ 

Given the estimates a and the estimate of labeling a(d) = 

(ai(d), a 2 (d), • • •,a|v|(d))'^ is determined by 

ai{d) = a.Tgm^^5a^xP^^\a\d) (*GV). (8) 

a 

This procedure to determine the estimate a{d) is approximately computed by using the 
LBP algorithm to 

We show numerical experiments by our proposed approach in Fig.2. In our numerical 
experiments, the test image d in Fig.2(a) is acquired from the Berkeley Segmentation 
Data Set 500 (BSDS500)^\ The size of the test image d is 321x481 and the sizes of d^^'^ 
are reduced to 20x30 for P = 8 and 10x20 for R = 10. The labeling configurations d 
obtained by means of our proposed algorithm based on our inverse RSRG transforma¬ 
tion for P = 8 and P = 10 are shown in Figs.2(b) and (c), respectively. In our proposed 
algorithm, and {m(^), C(^)|^gQ} for d^^'^ are estimated by using the conditional 
maximum entropy framework with the LBP such that they are obtained by subject¬ 
ing d^^^ to the algorithm presented in §3 of Ref.^fi The estimates for P = 8 and 
= 10 are 2.5288 and 2.5039, respectively. The estimates a = obtained by means 
of the inverse RSRG group transformations given by Eq.(7) for P = 8 and P = 10 
are 3.6765 and 3.6797, respectively. The labeling configuration in Fig.2(d) is computed 
directly by applying d to the algorithm in §3 of Ref.^\ The estimate P of a in Fig.2(d) is 
3.0301. The computational times corresponding to Figs.2(b)-(d) are 128(Sec), 78(Sec) 
and 1331 (Sec), respectively. Our numerical experiments were performed by using a per¬ 
sonal computer with an Intel(R) Gore(TM) i7-4600U CPU with a memory of 8GB. 
From the results of our experiment, we observe that the computational time can be 
reduced to less than one-tenth of that taken by conventional methods by application 
of our proposed algorithm based on the inverse RSRG transformation. We obtained 
results similar to the above mentioned case for other test images in the database set of 
BSDS500. 

In the present short note, we have presented a novel algorithm that involve the 
combination of the inverse RSRG transformation with the Bayesian image segmentation 
method proposed in Ref.^^. Our proposed method can reduce the computational time of 
the hyperparameter estimations significantly. We expect that our approach can also be 
applied to Bayesian image segmentations for three-dimensional computer vision, which 
remains one of significant problems. 


(a) (b) (c) (d) 



Fig. 2. Image segmentations in the case of 5 = 8 . (a) Test image d from the database of Ref.®^. 
(b)-(c) Labeling configurations by our proposed algorithm based on our inverse RSRG transformation 
for (b) i? = 8 and (c) R = 10. (c) Labeling configuration by Ref.^^ without using the inverse RSRG 
transformation. The results in (b)-(d) are shown with the color m(aj(d), d) at each pixel i for the test 
image d. 
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